Real space origin of temperature crossovers in supercooled liquids 



o 
o 

(N 



o 

£ 



13 

o 

> 
On 

\o 

O 

m 

O 

-i— > 
c3 



i 

c 

o 
o 



X 



Ludovic Berthier 1,2 and Juan P. Garrahan 1 

1 Theoretical Physics, University of Oxford, 1 Keble Road, Oxford, 0X1 3NP, UK 
Laboratoire des Verres, Universite Montpellier II, 34095 Montpellier, France 
(Dated: February 2, 2008) 

We show that the various crossovers between dynamical regimes observed in experiments and 
simulations of supercooled liquids can be explained in simple terms from the existence and statistical 
properties of dynamical heterogeneities. We confirm that dynamic heterogeneity is responsible for 
the slowing down of glass formers at temperatures well above the dynamic singularity T c predicted 
by mode coupling theory. Our results imply that activated processes govern the long-time dynamics 
even in the temperature regime where they are neglected by mode-coupling theory. We show that 
alternative interpretations based on topographic properties of the potential energy landscape are 
complicated and inefficient ways of describing simple physical features which are naturally accounted 
for within our approach. We show in particular that the reported links between mode coupling and 
landscape singularities do not exist. 
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I. INTRODUCTION 

The aim of this paper is to critically reconsider the 
physical origin of the onset of dynamical arrest and the 
associated crossovers between distinct dynamical regimes 
displayed by liquids supercooled through their melting 
temperature towards the glass transition 0, 0, 0, 01 • We 
do this by extending the real space theoretical framework 
based on dynamic facilitation of Refs. p| L4 ul 13 to the 
moderately supercooled regime corresponding to the re- 
gion where mode coupl ing th eory (MCT) supposedly 
applies, as reviewed in |lCtlll| . This novel approach takes 
directly into account the spatial aspects of the dynam- 
ics, in particular those related to dynamic heterogene- 
ity ^3 1 m contrast with many other theories 0, 0, 0] • 
Our analysis shows that the onset of slowing down can 
be understood in a simple physical way in terms of the 
dynamical properties of effective excitations, or defects, 
as a progressive crossover from a regime of fast dynam- 
ics dense in defect clusters, to one of slow heterogeneous 
dynamics dominated by isolated localized defects. Wc 
demonstrate that this real space picture explains the ob- 
served crossover temperatures, challenges the idea that 
these crossovers are related to changes in the topogra- 
phy of the energy surface or to MCT singularities, and 
is able to account for the apparent correlations observed 
between "landscape" and dynamical properties. 

The paper is organized as follows. In the rest of the 
Introduction we review the MCT and energy landscape 
points of view, discuss their problems and limitations, 
and describe the alternative real space perspective we 
will pursue. In Section[n]we develop the physical picture 
of the onset of slowing down and dynamical crossovers 
which emerges from our theoretical approach. In Section 
lllll we discuss its quantitative consequences and compare 
them to published numerical results. In Section IIVI we 
show how our approach also enables to derive the ob- 
served properties of the potential energy landscape of 
supercooled liquids. Finally, in Section we discuss our 



results and state our conclusions. 



A. MCT/landscape scenario... 

It is often assumed that the initial slowing down of 
the dynamics of supercooled liquids can be rationalized 
by MCT 0, 0, 0, US [ill. Numerical simulations are 
now able to investigate the first five decades in time of 
this slowing down [Tl|, so this is also the regime which 
has been studied in greatest microscopic detail. The 
degree of success of MCT is still a matter of debate. 
This is due to the fact that the central MCT predic- 
tion, a complete dynamical arrest at a temperature T c 
where the alpha-relaxation time diverges as a power law, 
t q (T) ~ (T — T c )~ 7 , is actually never observed, but 
a power law fit to the data apparently works on a re- 
stricted time window [ToL [lll | . The appearance of new 
mechanisms for relaxation, often termed "activated pro- 
cesses", but seldom described in any detail, is then in- 
voked to explain the discrepancy between observations 
and MCT predictions. In fact, activated processes are 
actually quantitatively denned, within MCT, by devia- 
tions between data and predictions 0, 0, Q|. It is 
believed that activated processes become relevant close 
to the dynamical singularity T c , their main effect being 
to prevent the predicted transition. 

From T c downwards it is assumed that the physics is 
dominated by activated processes, which determine also 
the canonical features of glass transition phenomena 1 1 ol | : 
non-exponential relaxation, strong and fragile liquid be- 
haviours, decoupling between transport coefficients, etc. 
It is sometimes said that the relevant physics for the glass 
transition sets in at T e , and is therefore out of reach of 
numerical simulations |16| . Crossovers into the activated 
dynamics reg ime are also reported to occur at tempera- 
tures T x [T3 or Tb depending on which aspect of 
the physics is considered. It is believed that all these 
temperatures are close enough to be taken as equivalent, 
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T C ^T X ^T B 0. 

The above scenario is apparently corroborated by the 
study of the statistical properties of the potential en- 
ergy landscape of model liquids 0, [23, 0, [S^. From 
the properties of the landscape two temperatures seem 
to emerge, T a and T c [23j. The onset of slowing down 
of the dynamics takes place at T , and coincides with 
the temperature below which the average energy of in- 
herent structures (IS), i.e., local minima of the poten- 
tial energy [2l|, e/s(T), starts to decrease markedly, see 
Fig. H This has been interpreted as the sign that the 
landscape starts to "influence" the dynamical behaviour 
|23j . At T c , it is further argued, a second change in the 
landscape properties takes place, which is indicated by 
several observations [H [H |H |H E3 |H El For 
example, the mean-square displacement from an equili- 
brated configuration to its corresponding inherent struc- 
ture, TV -1 J2i( r i ~ r i IS '') 2 j i s proportional to T below T c , 
as expected from pure vibrations in quadratic wells, but 
the temperature dependence changes above T c , revealing 
"anharmonicities" in the landscape [23l |24| , see Fig. ^ 

Another indication of a topological change in the en- 
ergy landscape was discussed in Ref. ISCl. in analogy with 
what happens in mean field models |3ll 132^ : the van- 
ishing as T approaches T c of the mean intensive num- 
ber of negative directions (intensive index) of stationary 
points of the potential e nerg y, n s (T) . Numerical simula- 
tions pi l3l l35l l3l I37L l3Sf found that n s {T) decreases 
with decreas ing T, and fits were performed to show that 
n s (T c ) = [33l 134. l35l|. The physical interpretation of 
this result is the apparent existence at T c of a "geometric 
transition" between a "saddle dominated regime" above 
T c and a "minima dominated regime" below T c |39|| . T c 
would then really coincide with the appearance of ac- 
tivated processes, described in a topographic language 
as "hopping" between minima of the landscape. Analo- 
gous findings had been previously reported from instan- 
taneous normal mode analysis of equilibrium configura- 
tions 40], also supporting a qual itative change in the 
landscape topology close to T c jU [43 • 

Figure summarizes this MCT /landscape scenario 
with published numerical data 0, [23[ 0] obtained for 
the standard supercooled liquid model of Ref. 0] . 



B. ... and its problems and contradictions 

At first sight Fig. ^ appears as convincing evidence in 
favour of the MCT/landscape interpretation of the dy- 
namics. On closer inspection, however, the above sce- 
nario is less robust. Several qualitative and quantitative 
observations do not fit into the picture presented above. 

(i) Activated dynamics above T c . The main idea behind 
the landscape approach [2(| is that vibrations and struc- 
tural relaxations take place on very different time scales, 
so that the system is "trapped" and vibrates in one min- 
imum before "hopping" to another minimum. This is 
indeed observed using the mapping from trajectories to 




FIG. 1: Onset of slowing down in the binary Lennard- Jones 
mixture of Ref. |T^. Three quantities are reported as a func- 
tion of temperature T. (i) Logarithm of the relaxation time, 
log r a (in arbitrary units extending over four decades in time) ; 
a MCT power law divergence, with T c — 0.435, was used to 
fit this data in Ref. [la. (ii) Energy of inherent structures, 
eis, taken from Ref. |43 |. which decreases markedly when 
the temperature decreases below T a = 1.0. (iii) Anharmonic 
part of Cartesian distance between configurations and their 
corresponding IS, (Ar) 2 = N' 1 £\(rj ~ r f S ') 2 - a T, which 
displays a qualitative change around T c , taken from Ref. |2^ |. 



IS [2l| in simulations of sufficiently small systems |44ll45| . 
We have also discussed theoretically this issue in a recent 
work 0. Given that numerical studies wereperformed 
much above T c , a crucial conclusion of Refs. 0,[4J,E3 is 
that "activated dynamics" is indeed present in this tem- 
perature regime. 

(ii) Heterogeneous dynamics above T c . It is now well- 
documented that the dynamics of supercooled liquids, 
even above the mode coupling temperature T c , is hetero- 
geneous in the sense that the local relaxation time has 
non-trivial spatial correlations [46j . This phenomenon 
is not very different from what happens experimentally 
close to the glass transition T g 12] . Besides, the de- 
coupling between transport coefficients, which is also in- 
terpreted in terms of dynamical heterogeneity > 
is observed in numerical simulations above T c [l3| . al- 
though the effect is quantitatively less pronounced than 
in experiments near T g |48|. 

(iii) Presence of saddles below T c . Despite claims 
based on numerical results that Tr marks a real change 
in the topology of the landscape [13, 13 HH 1 there are 
strong indications that this is at best only a crossover 
|36t l37| , and at worst a biased interpretation of numeri- 
cal data |38( . For instance, careful numerical studies have 
shown that the saddle index n s (T) remains positive even 
below T c [3g. Moreover, Ref. [33] argues convincingly 
that the data for n s (T) can be described by an Arrhe- 
nius law, n s (T) ~ exp(—E/T), with E an energy scale, 
which means that T c does not mark any particular change 
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in the saddle index. 



C. Alternative: real space physics and 
coarse-grained models 

The problems described above can be overcome 
through an alternative perspective on glass transition 
phenomena which puts the real space aspects of the dy- 
namics at its core. This is the approach developed in 
Refs. 0, IE 00 Interestingly several of its central con- 
cepts, like the relevance to the dynamics of localized ex- 
citations l49 l |50 | and the importance of effective kinetic 
constraints |^ Li3 > have been present in the literature 
for many years. Moreover, one of the original key ob- 
servations of dynamic heterogeneity in glass formers was 
mad e by Harrowell and coworkers |5^| in the models of 
Ref. See Ref. |3| for an exhaustive review. 

Our approach relies on only two basic observations, (i) 
At low temperature mobility within a supercooled liq- 
uid is sparse and very few particles are mobile. This is 
somewhat equivalent to the statement that particles are 
"caged" for long period of times, as reflected by a plateau 
in the mean-square displacement of individual particles, 
(ii) When a microscopic region of space is mobile it in- 
fluences the dynamics of neighbouring regions, enabling 
them to become mobile, and thus allowing mobility to 
propagate in the system. This is the concept of dynamic 
facilitation HD- The observation that very mobile 
particles in a supercooled liquid move along correlated 
"strings" |4(| is a confirmation of this fundamental idea. 

From these two concepts it is possible to build ef- 
fective microscopic models for glass formers by means 
of a coarse-graining procedure. This procedure can be 
schematically described as follows • Spatially, the par- 
ticles are coarse-grained over a length scale Sx of the 
order of the static correlation length given by the pair 
correlation function. This removes any static correla- 
tions between coarse grained regions of linear size 5x. 
Cells are then identified according to their mobility by 
performing a coarse-graining on a microscopic time scale 
St. In its simplest version cells are identified by a scalar 
"mobility field" , n(r, t) = 0, 1, the values 0/1 correspond- 
ing to an immobile/mobile cell at position r and time t. 
The next step is to replace continuous space by a lattice, 
n(r, t) — > rij(£). Mobile or excited cells carry a free energy 
cost, so when mobility is low it is reasonable to describe 
their static properties with a non-interacting Hamilto- 
nian [52j, 

N 

H = Y,n*, (1) 

i=l 

for a lattice of N sites. The link between mobility and 
potential energy [I| has also been observed in numerical 
simulations [46| . 

The coarse-graining procedure described above will 
generate local dynamical rules for the mobility field. The 



prominent feature of this dynamics will be dynamic facil- 
itation, which in its simplest version states that a cell at 
site i is allowed to move only if it has an excited nearest 
neighbour j^j, 

d c 

m = , m = 1, (2) 

Ci (1-c) 

where Cj = 1 — Y[<j a(1 — n j)i an d (j, i) indicates nearest 
neighbour, and c represents the average concentration of 
excited cells easily deduced from QJ: 

c(T) = {n i ) = {l + e 1 l T )-\ (3) 

Explicit examples where dynamic facilitation is generated 
under coarse-graining can be found in Ref. |55| . Clearly 
different models are defined simply by changing the ki- 
netic rules, e.g., the number or directionality of mobile 
neighbours required to move |54j. Also, a more complex 
mobility field may be required to account quantitatively 
for all glass transition features Q. 

Crucially, we will show that the physical mechanisms 
which explain the onset of slowing down and crossovers 
between different dynamical regimes in supercooled liq- 
uids are generic to this class of models. This means 
that we can use the simplest of them, the Fredrickson- 
Andersen (FA) model defined by Eqs. 11 121 in one spatial 
dimension (hereafter ID FA model) to make detailed pre- 
dictions and calculations. 



II. PHYSICAL PICTURE OF DYNAMIC 
CROSSOVERS 

In order to understand the physics captured by the 
coarse-grained facilitated models defined above, it is use- 
ful to look at trajectories, that is, space-time represen- 
tations of the dynamics 5]. We show in Fig. [21 three 
representative trajectories for the ID FA model, where 
mobile cells (defects) are black, and immobile ones are 
white. From the trajectories, the principal observation is 
the appearance at low temperatures of non-trivial spatio- 
temporal correlations, seen as spatially and temporally 
extended domains of immobile cells delimited by isolated 
defects In ID they look like "bubbles" 0, and tra- 
jectories are dense assemblies of these slow bubbles. This 
nanoscopic ordering in trajectory space is the cause of the 
phenomenon of dynamic heterogeneity observed experi- 
mentally and in simulations pj. Dynamic heterogeneity 
is the central aspect of the physics of supercooled liquids: 
it is naturally captured by our approach. 

The statistical mechanics of trajectories, rather than 
configurations, determines the dynamical behaviour. For 
example, due to the non-interacting Hamiltonian 
static correlations are trivial. However, when trajec- 
tories are considered, it is clear that cells become dy- 
namically correlated. In other words, these models nat- 
urally predict the existence of a dynamical correlation 




FIG. 2: Representative trajectories in the ID FA model. 
The vertical axis is space, the horizontal one time. The three 
trajectories are for L = 150 and t = 2000. Excited cells (or 
defects) are black, unexcited ones white. The top frame is 
for T = 2.5, in the high temperature regime where almost no 
isolated defects are present. The middle frame is for T — 1.0, 
the temperature regime where slow bubbles start to appear, 
seen here as large white domains. The bottom frame is for 
T = 0.5, where almost all defects are isolated. For T — 
0.5, the mean relaxation time is ~ 120, the mean dynamic 
correlation length ~ 9, but it is clear that times and lengths 
are broadly distributed. 



£{T) which grows when the dynamics slows down. This 
statement can be quantified 0, 0, 15(3. l57| by study- 
ing multi-point functions, for example, C(\i — j\,t) = 
(Pi(t)Pj(t)) - (P i (t)){P j (t)), where P,{t) is a dynamical 
correlator at site i (below we will consider the persistence 
of site i). The spatial decay of a function like C(\i — j\,t) 
defines unambiguously the dynamical correlation length 
£(T), as already discussed theoretically [ M B |5a| and 
measured in numerical simulations |46l l56t I57L l58t l59| . 
Furthermore, the joint distributions of time and length 
scales give rise to the canonical features of glass formers, 
like stretched relaxation, decoupling between transport 
coefficients, and (kinetic and thermodynamic) strong and 
fragile behaviours [EIS0- 

Let us take a closer look at the temperature evolution 
of the trajectories in Fig. [2 Starting from the very low 
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FIG. 3: Distribution of the logarithm of the persistence time 
of individual cells, 7r(lnr) for various temperatures. T = 1.0 
marks the appearance of a shoulder in the high temperature 
distribution. For T c < T — 0.6 < T , two "processes" coexist. 
Fast processes disappear close to T c = 0.3. 



temperatures where trajectories consist of a mixture of 
slow bubbles, the dynamics can be understood in terms of 
the opening and closing of bubbles, that is, the branching 
of an excitation line or the coalescence of two. As shown 
in |6(, these events are the "hopping between minima" 
described in Ref. [2(J. Therefore, we have a clear under- 
standing of "activated processes" and of their statistical 
properties 

As temperature is increased more and more defects are 
present. This has several consequences. First, the typ- 
ical spatial and temporal extension of bubbles reduces, 
that is, the system becomes faster and less heteroge- 
neous. Second, clusters of defects become more common. 
These objects are important because their dynamics is 
completely different from that of isolated defects. In a 
cluster, defects do not have to diffuse and create or an- 
nihilate other defects but can instantaneously relax in a 
much faster process. At high temperature, the dynamics 
is fast because almost no bubbles are present, and the dy- 
namics is governed by clusters of defects. Interestingly, 
at some intermediate temperature (middle frame in Fig. 
|2J) a coexistence between clusters and isolated defects is 
observed, so that the dynamics has a "mixed" character. 

Clusters of defects disappear much faster with decreas- 
ing temperature than the overall concentration of defects. 
The probability to have a cluster of k defects is indeed 
p(k) oc c fe , so that at low T we have p(l) ^> p(2) 3> ■ ■ • . 

The coexistence of fast and slow processes with differ- 
ent temperature behaviour has a direct influence on the 
distribution of local relaxation times, which we present 
in Fig. [21 for various temperatures. At high temperature, 
T> 1.0, where fast processes are dominant, the distribu- 
tion is exponential, with a mean which depends weakly 
on temperature (below we discuss in detail its tempera- 
ture dependence). Around T = 1.0, a shoulder develops 
in the large time tail of the distribution, corresponding 
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FIG. 4: Temperature regimes emerging from the discussion of 
section [21 T marks the onset of slow dynamics, the appear- 
ance of the isolated defects (bubbles, activated processes), and 
the growth of a dynamic correlation length. At T c , traces of 
the high-T physics (clusters) become negligible in the distri- 
butions of relaxation times. The crossover region T c < T < T 
has therefore a mixed character. 



cesses, while not dominant in number, still contribute to 
a significant fraction of the mean relaxation time (say 
again 90%) leads to the definition of T a : 



&t'-k{t')t' = 0.9 (r). 



(5) 



In Fig. 01 we summarize the physical picture that 
emerges from the considerations of this section. From the 
distributions of Fig. [3J we recognize that the dynamics 
becomes slow when the temperature is decreased below 
the onset temperature T Q . This distinguishes the trivial 
liquid and the slow glassy regimes. From the trajecto- 
ries of Fig. [21 we were also able to distinguish between 
fast, non-activated processes (clusters), and slow, acti- 
vated processes (isolated defects, bubbles). The former, 
typical of the liquid high-T physics, become negligible 
below T c , while the latter, typical of the low-T physics, 
appear at the onset temperature T D . As a consequence, 
the crossover region T c < T < T Q contains traces of both 
high-T and low-T physics, as observed in the time distri- 
butions of Fig. |21 



to the appearance of the bubbles in the trajectories of 
Fig. [3 This marks the increasing relevance of slow pro- 
cesses and the growth of the dynamic correlation length 
beyond the microscopic high temperature value. This 
temperature corresponds therefore to the onset temper- 
ature, T = 1.0. 

Decreasing further the temperature, we clearly see a 
regime of mixed dynamics. At T = 0.6, for instance, 
there are two peaks in the distribution, reflecting the 
coexistence of clusters and isolated defects. At this tem- 
perature, the time scale has already increased by several 
orders of magnitude, and the dynamic correlation length 
is about £(T = 0.6) « cT^T = 0.6) « 6. 

Finally, further decrease in temperature makes clusters 
of defects very rare and we are left only with the contri- 
bution of isolated defects. This low temperature distri- 
bution is the one discussed in Refs. which in turn 
implies the stretched exponential decay of dynamical cor- 
relators. The contribution of clusters becomes negligible 
beyond a second crossover temperature, here T c = 0.3. 
While this crossover temperature is not linked in any way 
to the mode coupling singularity T c , this choice of nota- 
tion will become clear shortly. 

From these distributions, it is possible to propose an 
empirical but quantitative determination of T D and T c . 
At each temperature T, the distribution is composed of 
fast processes, r < r*(T), and slow processes, r > t*(T), 
where r*(T) can be defined, e.g., as in Ref. [60|. Re- 
quiring that slow processes are a significant fraction (say 
90%) of the distribution leads to the definition of T c : 



(To) 



dr'7r(r') = 0.9, 



(4) 



where in an abuse of notation we also call 7t(t) the dis- 
tribution of persistence times. Requiring that slow pro- 



Ill. QUANTITATIVE CONSEQUENCES 

The physical picture we have presented for the onset 
of slowing down, based on the increasing relevance of a 
dynamically heterogeneous evolution of the system, leads 
to quantitative predictions which are in good agreement 
with previous numerical and experimental studies, as we 
discuss in this section. 



A. Dynamical correlators 

The basic dynamical quantities recorded in experi- 
ments and simulations of supercooled liquids are spa- 
tially averaged two-time functions. Simulations usually 
focus on the time domain and typically consider density- 
density correlation functions, while experimental results 
are often expressed in the frequency domain, measuring 
for instance dielectric susceptibilities. We will only con- 
sider systems in equilibrium so that the information con- 
tent of both kinds of measurements is equivalent. 

From the distributions of times, Fig. it is easy to 
derive dynamical correlators for the ID FA model con- 
sidered here. The spatially averaged persistence function 
reads, 



P(t) 



dr'Tr(r'). 



(6) 



The behaviour of P(t) as a function of time for various 
temperatures is shown in Fig. At very low tempera- 
tures, the persistence function is known exactly due to 
the diffusion properties of isolated defects, and one gets 



P(t) = exp 



-(T) 



(7) 
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FIG. 5: Persistence functions in the ID FA model for 
the same set of temperatures as in Fig. 01 Symbols arc 
numerical data and full lines are fits to the stretched ex- 
ponential form expected theoretically for low temperatures, 
P(t) = exp [~{t/T{T)f] , with /3 = 1/2. 



where t(T) is the relaxation time discussed in the next 
section. For the ID FA model, j3 = 1/2, but the stretch- 
ing exponent might be temperature dependent in more 
elaborated (fragile) models [EOSIj as is indeed observed 
in experiments [j], ■ 

We see from Fig.EJthat for T = 0.2 and 0.3 P(t) is well 
approximated by Eq. (0 on the whole time window. For 
higher temperatures, T > T c = 0.3, the mixed character 
of the correlators is evident from the fact that Eq. @ 
only describes the long time behaviour of the correla- 
tor, as expected. In this temperature regime, short times 
are best described by a simple exponential. In the high 
temperature regime, T > T Q = 1.0, relaxation is just ex- 
ponential for all times. We conclude that the appearance 
of isolated defects at T Q is reflected in the long time be- 
haviour of dynamical correlators. In the crossover region, 
T c < T < T Dl more and more of the decorrelation is due 
to isolated defects when T decreases. Below T c the entire 
decorrelation is due to these slow processes. 

A confirmation of the progressive domination of slow 
over fast processes described above can be found in the 
numerical results of Ref. [2|J. The similarity of Fig. 4 of 
Ref. |26| and our Fig. [S] is in fact quite striking. In par- 
ticular, Ref. j2|| calculated density-density correlations 
from both real configurations and their corresponding IS 
in a binary Lennard- Jones mixture. The latter, where 
thermal energies were removed in the quenching proce- 
dure, are the ones which have to be compared with Fig. 
|SJ since fast vibrations are also removed in our coarse- 
grained approach. 

An important conclusion is that the long-time decay 
of dynamical correlators, sometimes referred to as the 
alpha-relaxation, is due to the presence of isolated defects 
and therefore of heterogeneous dynamics, even in the 
T > T c regime. This means that activated dynamics, in 
the language of MCT, or hopping events in topographic 



FIG. 6: Temperature dependence of the relaxation time in 
the ID FA model, t(T). Open circles correspond to nu- 
merical data. Three fits are presented. The dashed line is 
the simple mean-field, Hartree like approximation, tmf ~ 
exp(l/T). The dotted-dashed line is the low-T exact be- 
haviour, r cx ~ exp(3/T). The full line is a power law MCT- 
like fit, t mct ~ (T - T c )~ 7 with 7 = 2.3 and T c = 0.3. The 
inset shows that the apparent power law behaviour is accept- 
able in a range of three decades in times. The main figure 
shows that low-T and high-T fits account for the whole tem- 
perature range. 

terms, are responsible for the alpha-relaxation, even at 
temperatures well above T c . This conclusion is unavoid- 
able in view of the numerical data of Refs. |26l B3 . l45| . 

B. Relaxation time 

The next natural quantity to consider, the relaxation 
time t(T), is readily obtained from the dynamical cor- 
relators discussed in the previous subsection. From the 
discussion of section [HJ we expect a crossover from high- 
T to low-T at the onset temperature T Q , see Eq. (JSJ. 
Our results for the ID FA model are presented in Fig. [ft] 
where t(T) is defined as the time where the persistence 
function has decayed to the value 1/e. 

The simplest mean-field approximation to the dynam- 
ics of the FA model consists in a Hartree like decoupling 
of spatial correlations, (fiiTij) — > rii{nj), in the dynam- 
ical equation for rij. This amounts to replacing the ac- 
tual neighbourhood of site i by an average neighbour- 
hood, and spins are always facilitated with an average 
rate equal to c. This gives a mean-field estimate of the 
relaxation time |6lj . 

tmf(T) » c- 1 ~ exp (I) . (8) 

Fig.EJshows that this simple approximation accounts for 
the dependence of the relaxation time at high tempera- 
tures, T > T . 

The exact result for the relaxation time of the ID FA 
model is obtained by realizing that isolated defects un- 



dergo diffusion with a temperature dependent diffusion 
constant, D(T) ks c ~ exp(— 1/T). The system relaxes 
when defects have diffused over a distance given by the 
mean separation between defects, c _1 , so that 



r cx (T) » D 



6X P I y 



(9) 



This mechanism relies on the notion of dynamic facilita- 
tion which implies that local fluctuations of the mobility 
determine the dynamics, and is essentially beyond reach 
of any mean- field type of approximation |53q . We see 
from Fig. that Eq. © accounts for the behaviour at 
low temperatures, T < T a , Figure HO also presents a fit to 
the data with an MCT power law form for the relaxation 
time p, 



tmct(T)«(T-T c )- 



(10) 



similar to the one obtained in Ref. 52] for the 2-spin 
facilitated, two dimensional version of the FA model. 

From Fig. we draw the following conclusions. The 
behaviour of the relaxation time changes from the high- 
T to low-T behaviour close to the onset temperature T Q . 
The combination of simple mean-field at high-T with the 
exact form at low-T allows to describe the temperature 
dependence of the relaxation time over the whole tem- 
perature range. However, given that t(T) smoothly in- 
terpolates between these two different functional forms, 
the MCT power law form, Eq. (JTUJ, appears to works 
reasonably well in a time window of about three decades 
(see inset in Fig. |BJ . This range of apparent power-law 
behaviour is in fact larger than the corresponding one in 
the canonical binary Lennard- Jones mixture of Ref. ^| , 
where extensive tests of MCT have been performed. Re- 
markably, we also find that the T c extracted from the 
power law fit to the relaxation time coincides well with 
the temperature where fast processes cease to contribute 
in a significant manner to the distribution of relaxation 
times, Fig. |3| This explains our choice of notation for the 
lower crossover temperature, T c . 

Following the standard MCT reading of the data 10, 
ITU Il3[ . we would erroneously conclude that activated 
processes only appear close to T c = 0.3, since these pro- 
cesses are often tautologically defined by the breakdown 
of the power law behaviour of the relaxation time 
Figs-Blllil and HO prove instead that activated dynamics 
starts to be relevant at T a , much above T c , dominating 
the alpha-relaxation of the correlators, and hence the re- 
laxation time of the system. The results of this section 
considerably weaken the possibility of the existence of 
a temperature regime in supercooled liquids where the 
relaxation time is correctly described by a power law be- 
haviour, ft follows that the standard determination of 
the location of the MCT "singularity" T r in experiments 
and simulations is physically unjustified [62j,|63j. In fact, 
the issue of the location of T c has been recently addressed 
in Ref. [64( i where it was found that for a variety of sys- 
tems, the temperature T c obtained from the actual MCT 
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FIG. 7: Onset of the slowing down in the ID FA model. 
We show three quantities as a function of temperature T. 
(i) Logarithm of the relaxation time, logr, see Fig. HJ (ii) 
energy of IS, e/s, which displays a qualitative change around 
T — 1.0; (iii) concentration d of cells moved in the descent 
from an equilibrium configuration to its IS, which displays a 
qualitative change around T c — 0.3. This figure should be 
compared with Fig. 



equations systematically coincides with the onset tem- 
perature T a discussed above. 



C. Crossover temperatures 

Let us now consider the quantities shown in Fig. 
as evidence in favour of the MCT/landscape picture in 
Lennard-Jones mixtures, from the perspective of dynam- 
ically facilitated models. In Fig. [7| we show for the ID 
FA model a plot analogous to Fig. ^ 

The first quantity presented in Fig. is the logarithm 
of the relaxation time, log r(T), as a function of tem- 
perature, which we discussed in detail in the previous 
subsection. 

The second quantity shown in Fig- His the average en- 
ergy of inherent structures, e/s. For the ID FA model 
it can be computed analytically b y s olving the zero- 



p y s ol 

temperature dynamics of the model [65|. 

e /s (T) = ce- c , 



(11) 



where the concentration of defects, c, is defined in Eq. 
At high temperature ejs changes very slowly. When T 
is reduced below T a the concentration of defects starts 
to decrease markedly, and e/s follows the same trend, as 
can be seen in Fig. [7] This change in behaviour at the 
onset temperature T Q is due to the appearance of isolated 
defects, and therefore of a heterogeneous dynamics, and 
not to any special change of the potential energy surface. 
This is a very different interpretation of the physics from 
that of Ref. 



The third quantity shown in Fig. is analogous to 
the distance between a configuration and its nearest IS 
(see Fig. 0). In the lattice models we are considering 
the natural quantity to compute is the concentration of 
sites which change during the descent towards the inher- 
ent structure, d(T). Since only excited sites can change 
during this procedure, we get, 

d{T) = c-e IS = c(l-e- c ). (12) 

Clearly, at low temperature d « c 2 ~ exp(— 2/T). This 
behaviour is physically natural. Contributions to d(T) 
come from clusters of defects, which are the only objects 
that can relax during the descent to an IS. Since the 
probability for a cluster of k defects, p(k), goes as p(k) f» 
c k , the main non trivial contribution to d(T) at low T 
comes from the smallest clusters, k = 2. These relax 
only one defect in the descent, so that d w p(2) « c 2 . 
Moreover, our defect interpretation is consistent with real 
space observations in simulations of silica j2^], where it 
was found that during the descent to the IS the major 
contributions to the distance comes from annihilation of 
localized topological defects of the amorphous structure. 

The similarity of Fig. [7| to Fig. is striking. The 
emerging physical picture is however completely differ- 
ent from the one of the MCT/landscape scenario. For 
example, while it may appear from the behaviour of 
d(t) above T c that this quantity extrapolates to zero 
when T — + T c (see Fig. 01, the exact temperature de- 
pendence of d(T) is purely Arrhenius. This means that 
T c has no particular importance [T c would not look spe- 
cial in a plot of d(T) versus l/T}. Even if one ac- 
cepts T c as delimiting two regimes with high and low 
concentrations of clusters of excitations, this apparent 
crossover is completely irrelevant as far as the long- 
time dynamics is concerned. These observations suggest 
that the crossover at T c from "landscape-influenced" to 
"landscape-dominated" of Ref. [2^] is not physically sig- 
nificant for the alpha-relaxation. 



IV. INTERPRETATION OF "LANDSCAPE" 
PROPERTIES 

In recent years, the potential energy landscape of su- 
percooled liquids has become an object of study per se 0. 
In Ref. Q , we have developed the idea that the main mo- 
tivation behind these works was the observation of the 
separation between fast vibrations and slow hopping pro- 
cesses if sufficiently small systems are considered. This 
apparently harmless statement on the system size, we ar- 
gued in Ref. |(| , results in fact from the central feature of 
the dynamics of supercooled liquids: "sufficiently small" 
really means "if the system size is of the order of the 
dynamical correlation length £(T)" [H,|25|]. However, in 
a purely topographic description of the physics based on 
the statistical properties of minima, the relevance of the 
dynamical correlation length is not obvious |66| . In that 
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FIG. 8: Top: zoom on the low-T trajectory of Fig. |U The 
vertical arrow indicates the closing of a bubble. Bottom left: 
expanded view of this event, showing two excitation lines 
meeting and coalescing. A cluster of three spins is needed 
for this process to occur. Bottom right: corresponding "re- 
action path". Before the event there are two isolated defects 
(energy = 2), a cluster of three defects (energy = 3) during 
the event, and one isolated defect after (energy = 1). 

sense, a topographic description of the glass transition 
misses a central aspect of the physics. 

We shall show below that using the very simple spatial 
approach described in previous sections, we can trivially 
derive the statistical properties of the landscape reported 
in recent years. This successful confrontation to such an 
amount of apparently non-trivial and detailed numerical 
results is again a strong indication of the validity of our 
approach. 

A. Real space description of "minima" and 
"saddles" 

Reference @ described in detail the connection be- 
tween the nanoscopic ordering in the trajectories of dy- 
namically facilitated models and the dynamics between 
IS, "metabasins" , or "traps" observed in numerical simu- 
lations or experiments of supercooled liquids. This same 
approach can be extended to account for the properties 
of "saddles" , i.e., configurations related to transitions be- 
tween IS. 

Figure|S|zooms on the lowest temperature trajectory of 
Fig. [21 The top panel of Fig.[S]shows diffusion of isolated 
defects. It also shows coalescence and branching events, 
i.e., closing and opening of bubbles. These two kinds of 
processes correspond to "hopping" events between dy- 
namical "traps" |fj . Lets consider in detail one of these 
events, for example the coalescence process enlarged in 
the bottom left panel of Fig. [8] Isolated defects diffuse 
by first facilitating one of their neighbours, for instance, 

10 -»• 11 ->• 01. (13) 



For two defects to coalesce the minimum number of exci- 
tations that have to be present when they merge is three: 



0.4 



101 ->■ 111 -> Oil ->■ 010. 



(14) 



In this sequence, the total number of defects is 2 at the 
beginning, 3 at the transition, and 1 at the end. This 
process is schematically described in the bottom right 
panel of Fig. [S] 

In topographic terms, an isolated defect corresponds 
to a local minimum of the energy, since such a configu- 
ration can only evolve by an energy increase, as in (|13fl . 
On the other hand, a cluster of three excitations corre- 
sponds locally to a saddle point, since it is the transition 
configuration between two minima, as in Q14f> and Fig. [S] 
Larger clusters thus correspond to higher order saddles, 
since the larger the cluster the larger the number of pos- 
sible moves into minima. The case k — 2 is particular: 
it is not a minimum since it can relax one defect to de- 
crease its energy, but it is not a saddle either since it 
does not correspond to a hopping event like that of Fig. 
IH1 Clusters with k — 2 are just ordinary points (i.e., not 
stationary points) of the landscape. The previous dis- 
cussion generalizes in a natural way to the whole class of 
dynamic facilitated systems. 



B. Absence of "geometric transition" 

The above identification between the relevant dynam- 
ical objects, isolated defects and clusters of defects, and 
"landscape properties" allows one to compute quantities 
such as the mean saddle index, n s (T), and the corre- 
sponding mean energy of stationary points, e s (T). 

The quantities n s and e s were estimated numerically in 
simulations of supercooled liquids [33L l34| . It was found 
that both functions decrease when T decreases, and ex- 
trapolations were performed that indicated n s (T c ) = 0. 
Also, plotting the dependence of n s on e s — ejs parame- 
terized by the temperature, a simple linear relation was 
obtained, n s oc (e s — ejs) [33. 134L l35|. 

In the case of the ID FA model, it is very simple to 
devise a procedure to go from an equilibrium configura- 
tion to the "nearest" stationary point. Isolated defects 
and clusters of defects with k > 3 are locally such points, 
so we only have to deal with k = 2 clusters. From these 
we can either reach a "minimum" k — 1 or a saddle with 
k = 3. We respectively assign the probabilities p s and 
(1 — p s ) to these two possibilities. We then have, 



,(T) = J2p(k)n s (k), 



k=l 



and 



e.(T) = X>(fc)e.(*), 



(15) 



(16) 




FIG. 9: Saddle index versus energy difference between saddle 
and minima computed analytically for the ID FA model for 
T £ [0, oo) and p s — 1/2. The obtained linear behaviour is a 
natural consequence of (i) dynamic facilitation, (ii) localized 
defects and (iii) dynamic heterogeneity. 



where p(k) = (1 — c) 2 c k is the probability to have a 
cluster of size k. From the discussion above we know 
that n s (l) = 0, e s (l) = 1, n.(2) = (1 - p.), e s (2) = 
p s + 3(1 — p s ), n s (k > 3) = e s (k > 3) = k. Putting all 
together we obtain, 



n s (T) = 3c 2 



(1 - Ps )(l - cY + c l--c 



and 



e s (T) = c+{l-2 P s)c 2 (l-cf 



At low T both quantities scale as n s 



and 



(17) 



(18) 



c. 



fc=i 



Three important conclusions can be drawn: 

(i) It is obvious from Eq. (HJ that n s (T) > for T > 0. 
This means that there is no "geometric transition" to a 
regime with vanishing saddle index. 

(ii) The saddle index has a temperature dependence 
which follows closely that of the distance d(T) discussed 
in the previous subsection. This is expected because they 
both receive their principal contribution, at low T, from 
clusters with k = 2. In other words, the main objects for 
the low temperature dynamics, the isolated defects, do 
not contribute to these quantities. Therefore, as for the 
distance d(T) in Fig. [3 the rapid decrease of n s (T) when 
the temperature decreases can easily be confused with a 
vanishing of the saddle index close to T c . 

(iii) The linear relation between n s and e s — e/s be- 
comes exact at low temperatures, see Eqs. I|11I17I18[1 . 
Again, the difference between e s and e/s comes from the 
clusters of defects which are relaxed during the descent to 
the inherent structure. In Fig. El wc show the behaviour 
of n s versus e s — e/s for the entire range T € [0, oo). Note 
that in Ref. |33| a linear behaviour between n s and e s was 
also reported. This is true at relatively high temperature, 
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given that above T a the energy of inherent structure is al- 
most constant while n s and e s change with temperature 
in the same way. 

These results are valid beyond the FA model which we 
have used to illustrate them. The inexistence of the geo- 
metric transition where the saddle index vanishes follows 
from the observation that the low temperature behaviour 
of n s is given by the smallest cluster of defects necessary 
to make a transition, in the sense described in the pre- 
vious section. Since these objects are spatially localized, 
their energy cost is 0(1), and they exist with non zero 
probability at finite temperature, T > 0. This argument 
is close in spirit to Stillinger's argument for the inexis- 
tence of an entropy crisis at the Kauzmann temperature 
involving point defects ■ Moreover, as discussed in the 
introduction, careful numerical simulations both confirm 
that n s (T < T c ) > [ME^ an d report an Arrhenius 
behaviour n s (T) |38j. in agreement with our results. 

The relation between saddle index and energy, 



oc (e s - eis), 



(19) 



which was first observed numerically [33|, yj, |3f|, is 
also a general result for dynamically facilitated systems. 
This relation contains two different pieces of information. 
First, it shows that the intensive saddle index is a number 
of 0(1). This is again a trivial consequence of the exis- 
tence of the dynamical correlation length £(T), so that 
a large sample can in fact be thought of an assembly of 
independent subsystems of linear size t(T) |(|. Second, 
and more interesting, is a connection between energy and 
saddle index, presented as a "general feature of the po- 
tential energy landscape of supercooled liquids" [3f| , for 
which no theoretical explanation was however available. 
This feature is in fact almost a tautology in the context 
of facilitated models: the more defects are present, the 
more available directions to move, the higher the energy 
above that of the IS. The fact that relation (|19f) holds 
in different model liquids is another confirmation that 
dynamical facilitation is a key generic feature of the dy- 
namics of supercooled liquids. 



C. Thermodynamics and "anharmonicities" 

Another common procedure of the landscape approach 
is to decompose configurations into vibrational and con- 
figurational components. Stillinger and Weber [22] sug- 
gested to perform this decomposition at the level of the 
partition function, 



Z(T) 



^(£ /s )exp 

E IS 



Eis 



F(T;Eis}\ 



T 



(20) 



where the sum is over energies of IS, Eis, their number 
is indicated by Q(Eis), and F(T, Eis) is the "basin free 
energy" which takes into account fluctuations within an 
IS due to vibrations and possible "anharmonicities" (i.e., 
all the rest). 



It is instructive to consider the calculation of the par- 
tition function in the case of the ID FA model using 
the Stillinger and Weber decomposition. The thermo- 
dynamics of the FA model is that of a non-interacting 
gas of binary excitations. This simple thermodynamics, 
however, can be obtained with any dynamics obeying de- 
tailed balance with respect to the Hamiltonian QJ, the 
actual FA dynamics defined by being just one possi- 
bility. In this sense, the Stillinger and Weber prescrip- 
tion for thermodynamics is an approximation for the way 
a thermodynamic quantity would be calculated using a 
particular choice of dynamics to sample the configuration 
space. Inherent structures, basin free energies, etc., have 
no thermodynamic meaning, they only have a dynamical 
meaning associated with a particular choice of dynamics. 

For the ID FA model we can evaluate the Stillinger 
and Weber partition function exactly. We have: 



N/2 



Zn{T) = 



Ers=0 



(7V - E IS )\ 
E IS \ {N-2E IS )\ 

1 



x exp ( -- [Ej S + F anh (T; E IS )} ) . (21) 



The first factor counts the number of configurations of en- 
ergy Ejs with only isolated defects in a system of N sites. 
Due to the coarse-grained nature of facilitated models the 
only contribution to the basin free energy comes from an- 
harmonicities. 

Performing the sum in Eq. (|21|l without the anhar- 
monic contribution gives, 



Z 



£T m (T) = exp 



N 
2T 



U N 



2 eXP l2^ 



(22) 



where Ujy(x) is the ZV-th Chebyshev polynomial of sec- 
ond kind. In the thermodynamic limit the above expres- 
sion simplifies to give the free energy in the "harmonic" 
approximation, 



/harm CO = nm ~^^Z N (T) 



T In 2 — T + vT+4e-V^)<23) 



In Fig. we compare the approximation (|23[1 to the 
exact expression for the free energy, 



/ex(T) = -Tin 1 



-l/T 



(24) 



As observed numerically in supercooled liquids [3j, both 
thermodynamic evaluations, Eqs. I|23|) and (|24|l . appar- 
ently coincide below T c when anharmonicities become 
negligible. From the previous sections we know that an- 
harmonicities are just a consequence of the existence of 
clusters of defects. At low temperature the difference 
between (|23f) and 124f> is therefore proportional to c 2 , 
reflecting the fact that: (i) anharmonicities do not disap- 
pear below T c , which again is no particular temperature 
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0.6 




FIG. 10: Comparison of the various expressions for the free 
energy (divided for convenience by —T) in the ID FA model: 
/ cx is the exact free energy I24H ; /harm is the purely harmonic 
evaluation 1231 which is a good approximation below T c = 0.3, 
and fi is the result obtained with the expression 1261 1 for 
anharmonicities which is good up to T = 1.0. 

in this context; (ii) anharmonicities are due to clusters 
of defects, k = 2 being the leading term at low tempera- 
tures. 

In the particular case of the FA model we can formulate 
an exact expression for the anharmonic free energy F^nh- 
It is easy to check that the choice 

Fanh(T; E IS ) = E IS {T) / CX (T) (25) 

in Eq. I|21|) yields the exact expression for the free en- 
ergy in the thermodynamic limit. This exact expression 
for the anharmonic part of the free energy is simple to 
understand. In an IS all defects are isolated. The con- 
tribution of clusters of defects is obtained from the prob- 
ability that an isolated defect in the IS was a cluster in 
the original configuration, which is given by Eq. I25J) . 

Interestingly, a numerical procedure to evaluate the 
anharmonic contributions can be proposed. First, the 
harmonic expression (|23[1 is evaluated. Then, the ap- 
proximate expression 

F a „h(T; Eis) ~ E IS {T) / harm (T) (26) 

can be used as an educated guess for the anharmonic 
contributions. This gives in turn a first order free energy, 
/i (T) . The improvement on the harmonic evaluation can 
be judged in Fig. 1101 where we see that fi{T) coincides 
with the exact free energy up to T ~ T Q . This evaluation 
can then be improved iteratively using F an h(T; Eig) ~ 
Ejs{T) /i(T) to get f2(T), and so on. It is easy to show 
that, in our particular case, limn^oc f n (T) — f cx (T). 

Although these results could lead to an improvement 
on present evaluations of anharmonic contributions in 
studies of the thermodynamics of supercooled liquids, 



they also show that topographic concepts are very far 
from the physical objects they pretend to describe. 



D. Failure of Adams-Gibbs relation 

We end this section with a remark on the Adam-Gibbs 
relation, which is an attempt to connect dynamical prop- 
erties to thermodynamic ones. The Adam-Gibbs for- 
mula relates the relaxation time of a glass former, r a , 
to the configurational entropy, S c [which would corre- 
spond to S c = (lnQ(Eis)} in the IS formalism]: r Q cx 
exp [l/(TS c )]. Apparently, this relation has been seen to 
hold both in numerical simulations and in experiments 
of various systems 3]. A careful look at the published 
data reveals however that the correlation between relax- 
ation time and entropy does not quantitatively satisfy 
the Adams-Gibbs relation. This important observation 
is often not clearly stated [3j. 

Our analysis shows indeed that thermodynamic prop- 
erties do not fully determine dynamical behaviours. 
Clearly, almost by definition, r a increases and S c de- 
creases as temperature is lowered, but that is where the 
connection ends. It is easy to check that the Adam- 
Gibbs formula fails completely when applied to dynamic 
facilitated systems. We find instead that time scales 
are broadly distributed, the distribution of times being 
the result of an integral over a distribution of length 
scales scale, p(£), imposed by thermodynamic equilib- 
rium, Eq. Crucially however, dynamics also enters 
the integral in the form of the conditional probability 
of time and length, p(t\£), which can be described, in 
a topographic language as containing informations on 
the relevant "barriers" , which have a priori no obvious 
link with the statistics of minima. As a consequence, 
thermodynamics alone cannot be used to predict the dy- 
namical behaviour. Again, we find in the literature an 
excellent numerical confirmation of this statement. In 
Refs 38, 60], using a purely topographic description of a 
supercooled liquid, it was shown that the diffusion con- 
stant could be computed by a combination of thermo- 
dynamic and dynamical quantities, well in line with the 
above discussion. 



V. CONCLUSIONS 

In this paper, we have developed a spatial description 
of the physics of the progressive slowing down of super- 
cooled liquids. The only ingredients in our method have 
been the notions of localized mobility excitations and fa- 
cilitated dynamics Our results were illustrated 
explicitly for the simplest case of the ID FA model, but 
are generic for this theoretical approach, and are in very 
good agreement with experimental and numerical obser- 
vations in supercooled liquids. 

The physical picture which emerges from our work 
is however markedly different from that of the 



12 



MCT/landscape scenario discussed in the Introduction. 

At high temperatures, T > T , the dynamics is fast 
and liquid-like, corresponding to the relaxation of large 
clusters of defects. Dynamic facilitation plays no major 
role, and a simple mean-field Hartree-like decoupling of 
the equations of motion yields predictions in good agree- 
ment with numerical results. 

When T < T a , the dynamics becomes heterogeneous, 
in the sense that local relaxation times are spatially cor- 
related in a non trivial way. This can be seen in the tra- 
jectories of Fig.|3as the appearance of slow bubbles 0,0 • 
The long-time dynamics of the system results from the 
wide joint distribution of length scales and time scales, 
and the relaxation becomes stretched. This dynamic het- 
erogeneity, which can be thought of as the activated dy- 
namics invoked, but never described, by MCT, deter- 
mines the alpha-relaxation and its temperature depen- 
dence for T < T a . Also, dynamic heterogeneity implies 
that decoupling of transport coefficients actually starts 
at T Dl as confirmed by the simulations. From a theo- 
retical point of view, local fluctuations of mobility cru- 
cially influence the dynamical behaviour. Any mean-field 
like approach, no matter how involved, is most probably 
doomed to fail. 

At T a not all trace of high-T physics (clusters of de- 
fects) disappears. The dynamics has a mixed character 
in the range T c < T < T a , as seen for example in the 
behaviour of dynamical correlators like in Fig. \5\ The 
temperature T c is just a crossover. It is the temperature 
below which isolated defects not only dominate the long 
time dynamics (as for T c < T < T Q ) but are also the most 
numerous dynamical objects, see Eq. flf. Clusters of de- 
fects, whose dynamics is homogeneous and non-activated, 
are responsible for the temperature dependence of sev- 
eral quantities, such as distance between configurations 
and IS, saddle index n s (T), and anharmonic contribu- 
tions to the free energy. In numerical simulations, T c 
has been interpreted as a key temperature, in accordance 
with MCT for which it represents a dynamical singular- 
ity. We have shown, however, that all of these quantities 
have a smooth temperature dependence, as has been re- 
cently observed numerically [3&|33,|38|]. This means that 
T c does not correspond to a transition or singular point, 
but is at most a crossover. Crucially, the objects which 
display a crossover close to T c are also irrelevant for the 
long-time dynamics, so that the inexistence of the singu- 
larity T c is anyway not a physically important issue for 
the alpha-relaxation. 

Below T c , isolated defects are the only remaining ob- 
jects and the dynamics is dominated by the nanoscopic 
demixing of slow and fast regions so that trajectories look 
like a dense mixture of slow bubbles, which in turn gives 
a natural theoretical interpretation of the canonical fea- 
tures of glass transition phenomena 0, IE • 

Our results, tog ether with some other recent stud- 
ies m h m hi m m m m m m, suggest that 

several essential features of the dynamics of supercooled 
liquids need to be recognized and we now list some of 



them. 

• The dynamics is heterogeneous and activated well 
above T c . 

• The dynamical slowing-down of supercooled liquids 
is due to the growth, below T Q , of a dynamic cor- 
relation length £(T), or more precisely, of a whole 
distribution of length and time scales. 

• The long-time dynamics, and therefore the relax- 
ation time r Q of the liquid is dominated by hetero- 
geneous "activated" dynamics below T a . 

• The MCT definition of activated processes as de- 
viations from the ideal theory is incorrect. It is 
unlikely that the power law behaviour predicted by 
MCT correctly describes the temperature depen- 
dence of r a . The practical definition of the tem- 
perature T c cannot be used. 

• No topological change of the potential energy land- 
scape takes place close to T c . Quantities such as 
the saddle index and anharmonicities do not van- 
ish close to T c and have a smooth temperature be- 
haviour. At best, they undergo a crossover from 
large to small which remains to be quantified. 

• Even if one accepts T c as a crossover temperature, 
as in Eq. quantities related to this crossover 
are unimportant for the long-time dynamics. 

• Knowledge of thermodynamic properties is not 
enough to predict dynamical behaviour, which ex- 
plains the quantitative failure of relations like the 
Adams-Gibbs formula. 

The approach we developed in this paper, which is an 
extension of previous efforts 01, El II] i is generic. It can 
be applied both to systems like Lennard-Jones liquids or 
to hard sphere systems. It gives a perspective on the 
physics of glass formers which is clearly distinct to, and 
in many respects much more natural than, that of MCT 
or topographic approaches. 

There are many important and interesting open ques- 
tions which need to be addressed from this new perspec- 
tive. This include, among others, understanding prop- 
erly the origin of mobility excitations, and the break- 
down of Stokes-Einstein-Debye relations and associated 
decouplings between transport coefficients. 

A general conclusion that can be drawn from this work 
and our previous ones is that, in many respects, glass 
transition phenomena is more standard than often as- 
sumed, in the sense that it is determined by the interplay 
between growing dynamic lengthscales and timescales. 
This is obviously reminiscent of critical phenomena |||, 
meaning that it should be possible to adapt renormaliza- 
tion group techniques to study the dynamics of the glass 
transition. 
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